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We study multi-qubit quantum channels that can be represented as a product of one-mode 
fermionic attenuation channels. An explicit formula for the classical capacity C\ and for the min- 
imum output entropy S m in of these channels is proposed. We compute S m in analytically for any 
number of qubits under assumption that the minimum is achieved on a Gaussian input. Apart from 
that, a simple numerical method for evaluating Smin is developed. The method is applicable to any 
channels that are sufficiently noisy. For fermionic product channels the proposed formula for Smin 
agrees with the numerical results with a precision about 10 -9 . 

PACS numbers: 



I. INTRODUCTION AND SUMMARY OF 
RESULTS 

Transmission of a classical information through a quan- 
tum communication channel is one of the basic problems 
studied in the quantum information theory. Suppose a 
sender (Alice) wants to send a string of n classical bits 
through a one-way noisy quantum channel <f> to a receiver 
(Bob). Suppose Alice is allowed to use the channel only 
m times. A number R is called an achievable rate if 
for any e, S > there exists a transmission protocol with 
n/m > R—S and an error probability smaller than e. The 
classical capacity of $ is defined as the minimum number 
C such that R < C for any achievable rate R. A closely 
related quantity is a one-shot capacity C\ which is de- 
fined analogously to C with a restriction that Alice uses 
only product states for signaling. It has been shown that 
C($) = linim^oo Ci(<i> lg " n )/m and conjectured that actu- 
ally C(<E>) = Ci ($), i.e., C\ is additive under tensor prod- 
uct of channels. The Holevo-Schumacher- Westmoreland 
theorem [lj, states that 

C 1 ($)=su PX (f) (1) 

£ 

where £ = {p a , p a € im(<I>)} is a probabilistic ensemble 
of states from the image of $ (a set of all output states) , 
and x(£) is the Holevo quantity 

\ a / a 

Important examples of quantum channels for which an 
explicit formula for the classical capacity is known are 
product of one-qubit unital channels 0] and products 
of bosonic attenuation channels Q. In the latter case 
a restriction on a power of the input signal has to be 
imposed to regularize the capacity. An explicit formula 
for the one-shot capacity is also conjectured for bosonic 
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attenuation channels combined with a classical Gaussian 
noise 

Quantum channels explored in this paper are direct 
fermionic analogues of bosonic product channels stud- 
ied in 4]. In contrast to their bosonic counterparts, 
fermionic modes are described by a finite-dimensional 
Hilbert space, so one does not need to regularize the ca- 
pacity. 

An algebra of observables of n fermionic modes for- 
mally coincides with the one of n qubits and can be con- 
veniently described by generators c\, c%, . . . , c.2n, 

c 2 j-i = cr{ (g) • • • <g) <S> of ® I ® • • • <8> J, 
c 2 j = of <8> • • • ® <8 crj ® / ® • • • <8 /. (3) 

Here cr" are the Pauli operators on the qubit j. This 
is the well-known Jordan- Wigner transformation that is 
often used to map a system of spins (qubits) into a system 
of fermions. The generators ci, . . . , C2n obey the Clifford 
algebra commutation rules 

^-p Cq &q Cp — 2 Spql , — Cp. 

An arbitrary linear operator acting on n qubits can be 
uniquely represented as a linear combination of 4™ mono- 
mials 

c(xj — c 1 c 2 • • • c 2n , 

where x = (x\, . . . , X2n) is a binary string of 2n bits. 
Define a quantum channel <!> by the following rules 

$(/) = I, 

$(c p ) = b p £p, p = 1, . . . ,2n, 
$(5(x)) - n ( 4 ) 

p : x p = l 

Here < b\, . . . , &2 rl < 1 are 2n real parameters specify- 
ing the channel (attenuation coefficients). A proof that 
the linear map defined in Eq. |0J is indeed a quantum 
channel (i.e. a trace preserving completely positive map) 
as well as some motivation of this definition will be given 
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in Section |HJ It suffices to mention now that <!> has a 
clear product structure. One can say that the system 
of n qubits is 'partitioned' into 2n separated subsystems 
(Majorana fermionic modes) that are described by the 
operators c%, . . . , ci n . Each subsystem is transmitted in- 
dependently through its own 'wire' described by a quan- 
tum channel c p — > b p c p . However this product structure 
should not be mixed with the tensor product of quantum 
channels. Indeed, we shall see in Section [H] that Krauss 
operators corresponding to different 'wires' do not com- 
mute. Accordingly, one can not expect that the capacity 
of the channel Eq. 10} will be additive (a capacity of an 
individual wire is not even a well-defined quantity). 

To compute the maximum at Eq. Q we shall use a cer- 
tain variational family of states, namely, fermionic Gaus- 
sian states [a, 13 ■ A Gaussian state p is completely speci- 
fied by its first and second moments Tr(pc p ), Tv(pc p c q ), 
see Section IlIII for a strict definition. In particular, the 
von Neumann entropy S(p) is a simple function of these 
moments. We will see that the product fermionic chan- 
nels map Gaussian states into Gaussian states. Define a 
Gaussian capacity as 

Cf($)= S up X (f 5 ), (5) 

where £ g = {p a , p a S 7ro($)} is an ensemble of Gaussian 
states (probabilities p a may be arbitrary though). This 
is the best transmission rate that can be achieved if Alice 
uses only product Gaussian states for signalling. Obvi- 
ously, Cf < C\. Our main result is an explicit formula 
for the Gaussian capacity. 

Theorem 1. Let $ be a fermionic product channel acting 
on n qubits with attenuation coefficients < b\ , . . . , bin < 
1. Denote bj the j -largest coefficient. Then 

where H(x) — — cclog(x) — (1— x) log(l— x) is the Shannon 
binary entropy. Moreover, Ci(<I>) = Cf (<&) ifn<2. 

A proof of the theorem as well as some additional re- 
sults that are valid for Gaussian input states is given in 
Section 

The theorem implies that the optimal signaling ensem- 
ble is Gaussian for n < 2. Can one achieve a higher trans- 
mission rate for n > 2 by using non-Gaussian ensembles? 
To get some intuition about it, we have tried to calculate 
C\ numerically for three and four qubits, see Section Ivl 
The numerical results strongly suggest that C\ = Cf, 
i.e., that we gain nothing from using non-Gaussian sig- 
naling states. It might seem rather surprising because 
mixed Gaussian states are known to be the maximally 
noisy states (in terms of their von Neumann entropy) for 
fixed first and second moments. On the other hand, one 
can easily see from Eq. (0J that the channel $ washes 
out correlations among large number of modes such that 



l-mode correlators acquire a factor b , where & is a typical 
value of b p . It means that two-mode correlations is the 
best place to keep the information that has to be sent 
through the channel. In this respect Gaussian states are 
very promising candidates, because they have only two- 
mode correlations. Therefore, there might be a subtle 
tradeoff between large entropy of Gaussian states and 
their special correlation structure that makes them opti- 
mal for signaling. 

In general, a computation of C\ (even numerical) is 
an extremely hard problem. Fortunately, the fermionic 
product channels possess a special symmetry property 
known as covariance, see 0, [tj. For covariant n-qubit 
channels the one-shot capacity is given by 

C 1 ($)=n-5 min ($), (7) 

where 

S mm {<P) =min5($(p)) (8) 
p 

is the minimum output entropy. Our derivation of the 
Gaussian capacity is actually a computation of the mini- 
mum output entropy for Gaussian input states. The anal- 
ogous problem has been already solved for bosonic Gaus- 
sian channels 0] and a conjecture has been made that 
the global minimum of the output entropy is achieved 
on a Gaussian input. Whether or not this conjecture is 
true in the fermionic world is an open question (which is 
equivalent to an equality C± = Cf). 

In Section E| we describe a simple algorithm that al- 
lows one to evaluate S m i n ($>) for many quantum chan- 
nels. Applicability of the algorithm is not restricted 
to the fermionic product channels. In general, the al- 
gorithm works well if <i> is not too close to an ideal 
channel. The main idea is to exploit positivity of the 
relative entropy and iteratively minimize S(&(\ip){ip\)). 
More specifically, we construct a sequence of pure states 
\ipo)i 1^2), ■ • •, such that \ipj+i) is the largest eigen- 
vector of an operator <5>*(log$(|-0 J )(-0 J |)). The initial 
state IV'i) is chosen randomly. One can easily show that 
S($(ipj + i)) < S($(ipj)), i.e., the corresponding sequence 
of output entropies is non-increasing. One should repeat 
the iterations until the output entropy keeps decreasing. 

In contrast to local search methods, such as the gradi- 
ent descent, the iterative minimization is able to escape 
from local minimums of the objective function. How- 
ever, one can not guarantee that the iterations converge 
to the global minimum, see Section for more detailed 
discussion. 

In Section [V]B we use the algorithm to find the min- 
imum output entropy of the fermionic product channels 
for three-qubit and four-qubit systems. The correspond- 
ing value of the capacity C\ matches the Gaussian capac- 
ity Cf with a precision about 10~ 9 which has the same 
order of magnitude as the numerical noise. It is therefore 
fair to make a conjecture that C\ = Cf in general. 

The rest of the paper is organized as follows. SectionlTTI 
contains all the necessary facts about fermionic product 



3 



channels. In Section IlIII we review some basic properties 
of Gaussian states. Computation of the minimal output 
entropy for Gaussian input is done in Section llVl Also in 
this section we show that the image of a fermionic prod- 
uct channel has a nicely ordered structure with respect 
to the majorization relation. The algorithm for compu- 
tation of the minimal output entropy and the numerical 
results obtained with its help are discussed in Section Ivl 



II. BASIC PROPERTIES OF THE FERMIONIC 
PRODUCT CHANNEL 

A. Krauss representation 

We start from constructing an explicit Krauss repre- 
sentation for the channel <!> defined in Eq. (0J. Consider 
a parity generator 

P=(-i) n ClC 2 ---C2n-iC2n- (9) 

In terms of qubits it looks as P = o\ ® ■ ■ ■ ® a^. One can 
easily check that 

c(x)P= (-l) |x| Pc(x), 

where |x| is the Hamming weight of the string x. Intro- 
duce Krauss operators 

K p = iPc p such that K p c q = (-l) Spq c q K p . 

A one-mode quantum channel 

%(P) = \{l + b p )p + ^{l-b p )K p pK p 

implements the transformation Eq. (0J with b q = 1 for 
all q ^ p. Note that <J> P is a trace preserving completely 
positive (TPCP) map iff \b p \ < 1. Since the one-mode 
channels commute, $ p o $ 9 = $ 9 o <£> p , we conclude that 

$=$ 1 O$ 2 0-'O$2„. (10) 

Accordingly, $ is a TPCP map for all \b p \ < 1. Ex- 
panding the product Eq. I|10|) yields the desired Krauss 
representation of $, 

*(/>) = £>(x)fi(x)pc(x)t, (11) 

X 

where 

2n 

pW = ^rIl( 1 + (- 1 ) |x|+ ^M- 

9=1 

The channels related by a transformation b p — ► —b p for 
some p are unitarily equivalent, so we can focus on non- 
negative values < b p < 1. 



B. Covariance 

The subsection summurazes some ideas proposed in 
the papers @, 0- Let S — {U\, U2, . . . , Ud} be a set of 
n-qubit unitary operators (which may or may not consti- 
tute a group) that completely randomize any quantum 
state: 

1 d 

d^2 u jp u j = I l T for an y state p- ( 12 ) 

3=1 

A quantum channel Q is said to be covariant if it com- 
mutes with all operators U € 5, i.e., 

®{UAU f ) = UO(A)U\ U eS, 

for any operator A. Covariance is a very useful prop- 
erty of a channel that allows one to reduce the problem 
of calculating the capacity to the problem of finding the 
minimum output entropy of a channel, see Eq. (JSJ. In- 
deed, let I 'J) be a state with the minimum output en- 
tropy. Consider an esemble of pure states {U\^)}ireS, 
where U is chosen randomly and uniformly from the 
set S. This ensemble has the maximally mixed average 
state since d" 1 Y? j=l Uj\^){^\U] = 2~" I. Therefore the 
Holevo quantity of this ensemble is n — >S'(©(|\E r ) = 
n — S m i n (Q). Obviously, this is the absolute maximum of 
the Holevo quantity for any ensemble of states in im(O), 
so that Ci(O) = n — S m i n (Q) for any covariant channel. 

The channel $ defined in Eq. Q is covariant. Indeed, 
let S be the set of d = 4™ operators c(x). One can easily 
check that 

X ^ 

It means that S possesses the randomizing property 
Eq. I)12H . On the other hand it is obvious that 

$(c(x)Ac(x) t ) = c(x)$(A)c(x) t 

for any operator A. Therefore <& is covariant channel, 
and its one-shot capacity can be found from Eq. Q. 



C. Physical motivation 

Under certain circumstances the fermionic product 
channels may adequately describe evolution of fermi sys- 
tems interacting with environment. As a toy model con- 
sider a fermi system that consists of a single fermionic 
mode (n = 1) which we describe by creation/annihilation 
operators a', a. Suppose the system interacts with 
an environment which also can be described by a sin- 
gle fermionic mode a\, a e . Suppose the interaction be- 
tween the system and the environment is just the hop- 
ping Hamiltonian Hi nt = g(a jf d e -\-a\a). We assume that 
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the initial state of the environment is maximaly mixed, 
p e = (1/2)7. The toy model channel is defined as 

= Tr e (Up®p e U^) , U = exp (~iH mt t). 

After straighforward calculations one gets 

$(/)=/, $(a) = Aa, $(a t a - aa f ) = A 2 (a t a - aa t ), 

where A = cos (gt). Introducing c generators ci = ai+a* , 
ic-2 = di — a\ we get: $(cj-) = Acj and <I>(ciC2) = A 2 C1C2. 
This is an example of the product fermionic channel 
Eq. J3J with n — 1 and 61 = 62 = A. More generally, 
if the environment possesses a superconducting order pa- 
rameter, the interaction Hamiltonian may include terms 
with particle-hole conversion, like in Andreev reflection 
of electrons. In this case one can tune the interaction to 
obtain an arbitrary string of coefficients b p . 

In fact, any unital channel that can be described by 
a quadratic interaction Hamiltonian between the system 
and the environment is unitarily equivalent to a channel 
from the family Eq. (JIJ, see Q for a proof. 

III. GAUSSIAN STATES OF FERMIONS 

In this section we describe a variational family of n- 
qubit states that will be used to minimize the output en- 
tropy of the channel Eq. (0J. These are Gaussian states. 
The simplest example of a Gaussian state is the Fock 
vacuum |0® n ) or any other product state diagonal in the 
standard basis. Denote 

_^ n _^ n 

p* = & n( j + vj) = ^ n^ 7 ix > ^-i^j). (13) 

i=i i=i 

An arbitrary Gaussian state can be converted into the 
standard form p\ by a unitary evolution with a Hamil- 
tonian quadratic in the operators c p . Here is a strict 
definition. 

Definition 1. A state p is Gaussian iff it can be repre- 
sented as 

p = Up x U\ U = exp (iJT 2 + Mi), (14) 

where H2 and Hi are Hermitian linear combinations of 
operators c p c q and c p respectively. 

One can easily check that any pure Gaussian state has 
a form U |0® n ), where U is as above. 

It will be more convenient to work with even Gaussian 
states which correspond to choosing Hi = in Defini- 
tion □ 

Definition 2. A state p of n qubits is even Gaussian iff 
it can be represented as 

p = Up x U\ U = exp(iH 2 ), (15) 

where H2 is a Hermitian linear combination of c p c q . 



One can easily check that any pure even Gaussian state 
has a form exp (1H2) |0®"). Our strategy will be to prove 
all statements first for even Gaussian states. Then we 
shall show how to establish a correspondence between 
Gaussian states of rt-qubit system and even Gaussian 
states of (n + l)-qubit system. 

Unitary operators U = exp (1H2) will be refered to as 
Bogolyubov transformations. Their conjugated action is 

2n 

U C p ^ ^ Hpq Cq, 

q=l 

where R is a rotation, RR T = I, det (R) = 1. Any rota- 
tion R S SO (2n) can be realized by a proper Bogolyubov 
transformation U. In the rest of the section we list some 
basic properties of Gaussian states. 



A. Wick's theorem 

As in the case of Gaussian probability distributions, 
an even Gaussian state p is completely characterized by 
its covariance matrix, 

1 

M pq = - 2 Tr IP ( a p a 9 " £ i £ p)} ■ 

All higher moments of p can be expressed in terms of 
M using Wick's theorem. Namely, for any even binary 
string x S {0, l} 2 ™, |x| = 21, one has 

Tr(pc(x))=*'P/(Af[x]), (16) 

where M[x] is a 21 x 21 submatrix of M obtained by 
selecting all matrix elements M pq for which x p = x q = 1 , 
and Pf stands for the Pfaffian of a matrix. For example, 

i~ 2 Tr(pcic 2 c 3 c 4 ) = M12M34 — -M13M24 + M14M23. 

If p is an even Gaussian state then all odd correlators 
vanish, 

Tr(pc(x)) = whenever |x| =21 + 1, 

which can be easily derived from the fact that p com- 
mutes with the parity generator Eq. (|5J|, Pp = pP . 

We shall use Wick's theorem to prove that the channel 
<& defined in Eq. Q maps the set of Gaussian states 
into itself. Also Wick's theorem will allow us to describe 
the action of $ as a simple transformation of covariance 
matrices. 



B. Admissible covariance matrices 

Recall that our goal is to use Gaussian states as varia- 
tional states to minimize the output entropy of the chan- 
nel Since a Gaussian state is completely characterized 
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by its covariance matrix, we have to find a set of admis- 
sible covariance matrices. Note that if M is chosen ar- 
bitrarily, the operator p defined by Wick's theorem may 
have negative eigenvalues, i.e., it might not describe a 
quantum state at all. 

Lemma 1. A real antisymmetric matrix M is a covari- 
ance matrix of an even Gaussian state iff 

M T M < I. (17) 

The corresponding state is pure iff M M = I. 

For a proof of the lemma see 0. Since M is a real 
antisymmetric matrix, its spectrum consists of n pairs of 
conjugated eigenvalues ±iXj, Xj > 0. We shall refer to 
the numbers Ai, . . . , A„ as singular values of M. The con- 
sistency condition Eq. l|17f) is equivalent to inequalities 
Xj < 1. The singular values Xj completely determine the 
spectrum of an even Gaussian state which is a product of 
n binary spectrums (1/2)(1 ± Xj), see Q. Accordingly, 
the von Neumann entropy S(p) = — Tr(plogp) is 

S( P ) = ±h( 1 -^-). (18) 

C. Reduction to even states 

A simple correspondence between the sets of Gaussian 
states and even Gaussian states can be established by 
adding one extra fermionic mode to the system (this idea 
was proposed by Knill ^2). Let p be an n-qubit Gaus- 
sian state. Consider a linear map £ that maps n-qubit 
states into (n + l)-qubit states according to 

£(p)=V p®(I/2)V\ V = exp(zjc 2n+1 ). (19) 

Here 1/2 stands for one maximally mixed qubit which is 
labeled by n + 1. Since V is a unitary operator, one has 

S(£(p)) = S(p) + 1 for any p. (20) 

Let us show that £{p) is an even Gaussian state for 
any Gaussian state p. Indeed, representing p in the form 
Eq. lfH)l one gets 

£(p) = (VU)p x ®(I/2) (VU)\ 

where U = exp {iH<i + iHi). Taking into account that V 
commutes with p\, we can rewrite it as 

£{p) = (VUV ] ) p\ <g> (1/2) (VUV^. (21) 

Recall that H2 is a linear combination of c p c q with 1 < 
p,q < 2n, while Hi is a linear combination of c\, . . . , Qu- 
it follows that H2 commutes with V. On the other hand, 
V does not commute with Hi. Taking into account that 

V CpV^ = ic p C2„+i for any p = 1, . . . , 2n, 



one can easily check that V H{V^ is a linear combina- 
tion of operators c p c.2 n +i, P — 1> . . . , 2n. It means that 
V(H2 + H\)V' is a linear combination of quadratic op- 
erators c p c q only, where 1 < p, q, < 2n + 1. Comparing 
Eq. I|21l) and Definition [5] we conclude that £(p) is an 
even Gaussian state. 

Conversly, if 77 = £(p) is an even Gaussian state then p 
is a Gaussian state. Indeed, by definition of £, the gener- 
ator C2n+2 commutes with rj. It means that a covariance 
matrix M of r\ has all zeroes in the column 2n + 2 and 
the row 2n + 2, i.e., M p ^ n +2 = M 2n +2, P = 0. Therefore 
it can be represented as 

( N 0\ 

M = RM R T , M = 0, 
\ 0/ 

where R is a Bogolyubov transformation that involves 
only the generators £1, . . . , C2 n +i- Accordingly, r\ can be 
represented as 

r ! = Ui l0 ®(I/2)U\ U = exp(iH 2 ), 

where r/o is an even Gaussian state of n qubits with the 
covariance matrix N, while H2 is a Hermitian linear com- 
bination of c p c q with 1 < p,q < 2n + 1. Repeating the 
same arguments as above we get 

p® (1/2) = Wtjo ® (J/2) W\ W = V^UV. 

Conjugating H2 by V we do not affect terms c p c q that 
do not contain C2 n +\- On the other hand, a term like 
c P C2n+i is converted into ic p . Thus the generator C2 n +i 
does not enter into W, and we get p = W r/o W . It follows 
that p is a Gaussian state. 

D. Action of product channels on Gaussian states 

Let p be an even Gaussian state with a covariance ma- 
trix M and $ be a product fermionic channel with at- 
tenuation coefficients bi, ... , &2n- Moments of the state p 
are related to moments of the state $(p) by the following 
identity 

Tr(a>(p)c(x)) = j J] b p ) Tr(pc(x)). 

Comparing it with Wick's theorem Eq. Ijlfijl we conclude 
that <&(p) is an even Gaussian state with a covariance 
matrix 

<P(M) = BMB T , (22) 
where B is a diagonal matrix with entries b±, ... , &2n > 
B = diag(6i,6 2 , ■ ■ ■ Mn-iMn)- 

Lemma 2. Fermionic product channels map Gaussian 
states into Gaussian states. 
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Proof. Let $ be a channel as above. Define a fermionic 
product channel $ acting on n + 1 qubits such that 



$(cp) = fopC p , p = 1, . 

*(c 2 n+l) = C 2n+ i, 
*(C2„+2) = 0. 



.,2n, 



(23) 



One can easily check that £ commutes with $ in the 
following sense: 



£ o $ = $ o £ . 



(24) 



Let /j be a Gaussian state. Then £(p) is an even Gaussian 
state, and thus ($ o £)(p) is an even Gaussian state. It 
follows from Eq. (|24|) that (£ o $)(/?)) is also an even 
Gaussian state. As we proved above, it implies that $(p) 
itself is a Gaussian state. □ 



E. Three-qubit states 

Suppose we restrict our attention only to even states, 
i.e., those satisfying 



where Xj are real numbers and R G SO(6) is a rota- 
tion. Consider a Bogolyubov transformation U such that 
UdpW = J2q=i R pq Cq- Then 

p' = U pU ] = -(I + P)(I - i\i cic 2 - iX 2 c 3 c 4 - iX 3 c 5 c 6 ). 



Here we used the fact that UPW = det (R)P = P. The 
state p' is diagonal in the standard basis. On the other 
hand p 1 is a pure state. Therefore p 1 is a vector of the 
standard basis with even number of ones. One can easily 
check that any such vector is a Gaussian state. It follows 
that p is a Gaussian state as well. □ 

The minimum number of qubits supporting even pure 
non-Gaussian states is n — 4. For example, consider a 
state 

\6) = cos (0) |0, 0, 0, 0) + sin (0) |1, 1, 1, 1). 

One can easily check that the covariance matrix of |0) has 
singular values smaller than one whenever sin (20) =/= 0. It 
follows from Lemma^that such states are not Gaussian. 



IV. COMPUTATION OF THE GAUSSIAN 
CAPACITY 



where P is the parity generator, see Eq. @. Obvi- 
ously, for one qubit, n = 1, the only even state is |0). 
This is a Gaussian state. If one takes two-qubits, even 
states constitute a two-dimensional subspace with a basis 
|0, 0) and |1, 1). Let us consider operators A = (— i)c\c 2l 
B = (-i)c 2 c 3 , and C = (-i)AB = (-i)c 3 ci. They all 
commute with P and obey the same commutation rules 
as the Pauli operators a x , a v , a z respectivly. Therefore 
any unitary operator on the even subspace can be rep- 
resented as U = exp (iaA + i(3B + ijC) for some real 
numbers a, (3, 7. Accordingly, any even state has a form 
U 1 0, 0) . It follows from Definition [21 that any even two- 
qubit state is Gaussian. Surprisingly enough, the same 
is true for three qubits as well. 

Lemma 3. Any even pure state of three qubits is Gaus- 
sian. 

Proof. Let p = |^ r )(^| be an even pure Gaussian state of 
three qubits. Taking into account that all odd moments 
of p vanish and that P\^>) = \^>), we can express p as 



where the sum runs over all 1 < p, q < 6 and M is a real 
antisymmetric matrix. It is a well-known fact from linear 
algebra that such a matrix M can be represented as 



M = R 



Xiuj 
X 2 lu 
X 3 uj 



R 1 



-1 

1 



Throught this section $ is a fermionic product channel 
with attenuation coefficients b\ , . . . , 6 2n ordered in non- 
increasing way, i.e., 

h > b 2 > ■ ■ ■ > b 2n - 

A computation of the Gaussian capacity will proceed in 
three steps. Firstly we shall compute the minimum out- 
put entropy achievable on even Gaussian input states. 
Then we shall generalize it to arbitrary Gaussian input 
using the trick with addition of a qubit. Finally we shall 
verify that the optimal signaling ensemble consists of 
Gaussian states. 



A. Even Gaussian input 

Our intermediate goal is to find the minimal value of 
the output von Neumann entropy S(3>(p)) provided that 
the input p is an even Gaussian state. Denote this mini- 
mal value S min , e (<&). 

Let M be the covariance matrix of p. Since the min- 
imum output entropy is achieved on a pure input state, 
we can assume that 



M T M = L 



(25) 



see Lemma ^ As we have shown in Section ITTTl the out- 
put state <&(p) is an even Gaussian state with a covariance 
matrix 



$(M) =BMB T , B = diag(6i 



(26) 
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Let Ai, . . . , A„ be singular values of the matrix $(M) (re- 
call that its eigenvalues are Taking into account 
Eq. ||IEJ we get 



Smin,e = min ^ H 

3=1 



1 + A, 



where the minimum is taken over all strings of singular 
values Ai, . .., A n that are consistent with Eqs. Q25I26|) . 
We claim that 



3=1 



1 + hj-ifoj 



(27) 



The proof will be based on convexity arguments and two 
basic facts from the majorization theory, see [l^. Recall 
that if y = (yi, . . . ,y n ) and z = (z\, . . . ,z n ) are strings of 
real numbers, the majorization relation y -< z is equiva- 
lent to inequalities 



k k 
3=1 3=1 



(28) 



with an equality for k = n. Here yj and zj are the j- 
largest elements of the strings y and z respectively. 
Fact 1: If j(x) is a concave function of a real variable x 
and (yi, . ..,y„) -< {z\, ...,z n ) then 



3=1 3 = 1 

Fact 2: (Horn-Visser-Zaanen theorem) 

For any square matrix A define Di(A) as a product of I 
largest singular values of A. Then for any mxm matrices 
A and B 

Di(AB) < Di{A) Di(B) for all l = l,...,m. (29) 

Note that D m (A) = | det (A) | , so Eq. becomes an 
equality for I = m. 

Now let us proceed to a proof of Eq. I|27|l . Let Aj be the 
sequence of singular values Xj ordered in non-increasing 
way. Taking into account that 



Di(B) = Di{B T ) = \{ bj, Di{M) = 1, 

3=1 

we can rewrite the Horn-Visser-Zaanen inequality 
D 2k (BMB T ) < D 2k {B)D 2k {M) D 2k (B T ) 



k 2k 
3=1 P=l 



(30) 



with equality for k = n. Let us introduce auxiliary vari- 
ables yj — log (Aj) and Zj — log (6 2 j-i&2j) taking values 



on the interval (— oo, 0]. Obviously, y\ > . . . > y n and 
Z\ > . . . > z n . The inequalities Eq. I|30|) become equiva- 
lent to the majorization inequalities Eq. (|28|l . so that 



(yi, ...,y n ) < (zi, . . .,z n ). 
The von Neumann entropy of the output state is 

l-e x 



3=1 

To prove that the function f(x) is concave let us represent 
it as f(x) = H(g(x)), where g{x) = (1/2) (1 - e a ). Obvi- 
ously, 5(2;) is concave and maps the interval (— 00, 0] into 
the interval [0,1/2). The function H(x) is concave and 
monotone increasing on the interval [0,1/2]. Therefore 
H(g(x)) is concave on the interval (— 00, 0]. Applying 
Fact 1 we get 



3=1 3=1 

3=1 



(31) 



This value of the output entropy can be achieved if the 
input p is the Fock vacuum, i.e., 

1 - 

P=^II( 7 -^3-ic 2 3) = |O ")(O 8 "|. 
3=1 

Indeed, the corresponding output state is 
1 " 

®(P) = ^Yli 1 - ib 2 3 -lb2jC 2j -lC 2j ) 

3=1 
n 

= ^n^ + ^-i^crl). (32) 

(33) 



3=1 



The von Neumann entropy of <fr(p) matches the lower 
bound Eq. (|3T|) . 



B. Proof of Theorem [T] 

Now we are ready to prove Theorem^ Firstly we shall 
find the minimum output entropy achievable on a Gaus- 
sian input (not necesserily even). Denote this quantity 
S m in,g(&)- We claim that 

w*)=^( i ^)+E^( 1± ^ ±1 

(34) 

Indeed, let G n and be sets of n-qubit Gaussian 
states and even Gaussian states respectively, such that 



S 



min, g = min S($(p)). 
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We shall make use of the linear map £ defined in Eq. i|19|) 
that adds one qubit to the system. Taking into account 
Eq. JU we get 

Smin.g = min S((£ o $)(p)) - 1. 

Applying the commutation rule Eq. I|24|) we arrive to 

Smin, g = min S(($ o £){p)) - 1. 

peG„ 

As we have shown in Sect ion ITTTl C. £ (p) is an even Gaus- 
sian state, so that 

> min S(<t(j])) — 1. 

Since we already now how to compute the minimum over 
even Gaussian states, see Eq. (j57J), one gets 

S m in,q ^ / , H I — 1. 

U v 2 / 

Here bi, . . . , b 2n+2 are the attenuation coefficients for the 
map $, see Eq. (I23[l . Obviosuly, the largest of them is 
b{ = i>2n+i = 1, while the smallest one is S 2jl+2 = b 2n+2 — 
0. Therefore we arrive to 

S min , g > H (i±^l) + £ H + . (35) 

This lower bound is achieved on a state 

P* = K^i 1 + - «C 2 C 3 ) •••(/- «C 2 „_ 2 C 2n -l). (36) 

2" 

We have to verify that p* is indeed a Gaussian state in 
the sense of Definition ^ Indeed, one can easily check 
that 

p* = (v^io^ko^k^v^, 

where V 2 is a Bogolyubov transformation that shifts the 
operators c p cyclically according to 

(ci, c 2 , c 3 , . . . , c 2n ) — * (-c 2 „, ci, c 2 , . . . , C 2n -l), 

while Vi = exp (— ij c 2n ). The operator VI is chosen such 
that V\{c2nC\)Vi = —ic\. One remains to observe that 

^ n—l 

$ (/>*) = ^( 7 + b l 6 l) JJi 1 ib 23 b 2j+1 c 2j c 2j+1 ). 

i=i 

All factors in this product can be diagonalized simultan- 
iously, so the spectrum of $(p*) is a product of n binary 
spectrums (1/2)(1±6 X ) and (l/2)(l±& 2i 6 2j+1 ). The cor- 
responding entropy S(&(p*)) matches the lower bound 
Eq. 



It is obvious that Cf < n — S m ; m .g. The optimal sig- 
naling ensemble that achieves this upper bound consists 
of 2 n Gaussian states 

-^(1 ± ci)(I ± ic 2 c 3 ) ■■■(! ± ic 2 „_ 2 c 2n _i). 

These states constitute an orthonormal basis and each of 
them yields the output entropy S m in,g- 

To complete the proof of Theorem ^ one remains to 
show that C\ = Cf for n — 2, or, equivalently, that 
Smin — S m i n>g for two-qubit channels. Denote D n and 
sets of all n-qubits states and even states respectively. 
Repeating the same arguments as above one gets 

min S($(p)) > min S(^(rj)) - 1. 

pED 2 V&Dl 

Obviously, the minimum in the righthand side is achieved 
on a pure state r\. But we have already shown in Sec- 
tion [^^D that any pure even three-qubit state is Gaus- 
sian. Therefore one can substitute the minimum over D\ 
by a minimum over G\ and proceed as above. We have 
proved Theorem 

C. Minimal entropy and majorization 

In the previous subsection we have found a Gaussian 
state p*, see Eq. lj3o]l. such that S{<&{p*)) < S($(p)) 
for any other Gaussian state p. Since the von Neumann 
entropy measures amount of randomness contained in a 
state, we can expect that §(p*) is the least randomized 
state in the image of $ (as far as Gaussian inputs are 
concerned) . Here we shall put this statement into a more 
strict form and prove the following interesting fact. 
Proposition 1: 

4>(p) -< $(p*) for any Gaussian p. (37) 

Here the majorization relation between the density 
operators p < r\ means that the spectrum of p is ma- 
jorized by the spectrum of r\. The majorization relation 
p -< r] is the quantative version of the statement "p is 
more randomized than 77" . It is well known that p -< r\ 
iff P — ^2 a PaU a T]U^ for some unitaries U a and some 
probability distribution p a . In particular p -< ?; implies 
S a (p) > S a (ri), where S a (p) = T ^log(Trp a ) is the 
Renyi entropy. 

Proposition 1 is a simple consequence of the inequality 
Eq. i|30|) and the following technical lemma. 

Lemma 4. Let a\, . . . , a n and /3i, . . . , j3 n be real num- 
bers satisfying inequalities 

1 >«!>...> q„ > 0, 1 > /?i >...>/?„> 0, 

Let x = (x±, . . . ,x n ) and y — (yi, . . . ,y n ) be random n- 
bit variables with probability distributions 

1 ™ 

p w= ^ri( i +(- i ) a:j ^)' 



9 



and 



1 - 



If for all k = 1, 



3=1 
, n one has 

k k 
3=1 3=1 

wii/i an equality for k = n then P is majorized by Q, 

Let us first show how to derive Eq. (|37|l from the 
lemma. 

Proof of the proposition. It will be convenient to use a 
notation P n (a\, . . . ,a n ) for a product of n binary distri- 
butions (1/2)(1 ± ctj). One suffices to prove Eq. (|37l) for 
pure states p. As we have shown in the subsection A, 
<£>(/9) is an even Gaussian state whose covariance ma- 
trix has singular values A j , . . . , obeying the inequal- 
ity Eq. (|3L)[1 . As was mentioned in Section IIIII B. the 
spectrum of $(p) is a product of n binary spectrums 
(1/2) (1 ± Aj). Applying Lemma E| to otj = Aj and 
[3j = b 2 j-\b 2 j we conclude that 



$(p) -< P n (bib 2 ,...,b2n-lb2n)- 



(38) 



Now we use the correspondence between Gaussian n- 
qubit states and even Gaussian (n + l)-qubit states. Let 
p be an arbitrary Gaussian state. Consider a linear map 
£ defined in Eq. lfTO|l and a channel $ defined in Eq. (j2"3*|l . 
Recall that £ o $ = <f> o £. A state 

(| =(£o$)(p) = (lo%) 

is even Gaussian. Applying Eq. lj3l3|) to the channel $ we 
get 

Spec(rj) -< P n+ i(b 1 ,b 2 b 3 , . . . , 6 2 n-2&2n-i, 0). (39) 

Here Spec(ji) is the string of eigenvalues of r/. On the 
other hand, the spectrum of 7/ and the spectrum of §(p) 
are related by 

Spec(n) = Spec($(p)) X {-, -}, 

since £ just adds one maximally mixed qubit to $(p) and 
then applies a unitary operator that does not change the 
spectrum. Taking into account Eq. (|39f) we get 

Spec(<Z>(p)) x{i l}< 

P n+ l(b 1 ,b 2 b 3 , . . . , 6 2 n-2&2rc-l,0) 
= P„(&l,&2&3,...,fr2n-2&2n-x) X {\,\) ■ (40) 



This majorization relation involves 2 n+1 inequalities. If 
one picks out only 2" even inequalities (i.e. those includ- 
ing a sum of two, four, six, e.t.c. largest eigenvalues) one 
gets exactly the majorization relation 



Spec($(p)) -< P n {bi,b 2 b 3 , b 2n - 2 b 2n -i). 
It is equivalent to the statement of the proposition. 



□ 



Proof of Lemma We shall first prove the lemma for 
n = 2 by brute force method. Then we shall use standard 
majorization theorems to extend the proof to all n > 2. 

If n = 2 then the ordered probability distributions P = 
{pi,..-, pi} and Q = {qi,...,q{} are 



lA = 3(l + ai)(l + a 2 ), 
p l 2 = i(l + ai)(l-a 2 ), 
Pi = 3 (l-a 1 )(l+a 2 ), 
|(1 -<*!)(! -a 2 ), 



. = i(1 

4 = la 



/3i)(l + /3 2 ), 

q l 3 = |(l-/?i)(l + ft), 
qi = 1(1-/30(1-/32). 

The condition of the lemma for n = 2 looks as 



Pi 



c*i < (3i, ai« 2 = /3i/?2- 

Our goal is to prove that P -< Q. After a simple algebra 
the first majorization inequality < can be written 
as 

C C 
ai + — <Pi + —, C = a 1 a 2 =Pi0 2 - (41) 
Oil Pi 

It is satisfied because a function f(t) = t + C/t is 
monotone increasing on the interval t > \/~C and be- 
cause tt\,fi\ > \[C. The second majorization inequality 
Pi + P 2 < 9i + <7 2 i s equivalent to «i < f3\. The third 

majorization inequality +P3 < <7i+<72+<73 reduces 

to Eq. (j4T|) that we have already proved. Thus P ~< Q. 

To prove the lemma for arbitrary n > 2 we shall need 
two more basic facts from the majorization theory 12} : 

Definition 3. Let y = (yi, . . . , y n ) and z = (21, ... , z„) 
&e strings of real numbers. A string y is a T -transform 
of string z iff there exist a real number < A < 1 and 
integers 1 < a < b < n such that 

Vc = z c for all a, 6, 



Va 




A 1 - A 




Za 


Vb _ 




1 - A A 




Zb 



Fact 3: y -< z iff z can be converted to y by a sequence 
of T-transforms. 

Fact 4: y -< z iff ?/j = Pjk z k for some doubly 

stochastic matrix P. 

Consider auxiliary variables yj = log(<x/) and Zj = 
log(/3j). The condition of the lemma is equivalent to 
y -< z. Keeping in mind Fact 3, one suffices to prove 
that P -< Q whenever a string y = (7/1, . . . , y n ) is a T- 
transform of a string z = [z\, . . . , z n ). Let a, 6, A be as in 
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Definition Since a T-transform is a doubly-stochastic 
matrix, one has 

(2/o,2/6) -< (Zo,Zb). 

In terms of the original variables a and j3 it means that 

max(a a ,a b ) < max (f3 a , /3 b ), a Q a b = /3 a /3 b . 

Applying the n = 2 proof given above to the marginal 
probability distributions P ab and Q a b of the bits a, b we 
conclude that P a b -< Qab- Besides, for any bit c ^ s,f) 
the marginal distributions P c and Q c coincide. Fact 3 
implies that Q a b can be converted to P ab by a doubly 
stochastic matrix. Therefore Q can be converted to P by 
a doubly stochastic matrix and thus P < Q. □ 



V. NUMERICAL SIMULATIONS 

As we have seen before, the smallest number of qubits 
for which the Gaussian analysis fails is n = 3. Therefore, 
we shall carry out numerical simulations for fermionic 
product channels acting on three and four qubits. Our 
goal is to verify whether the minimum of the output en- 
tropy is achieved on a Gaussian input state. Below we 
describe an algorithm that in many cases allows one to 
find the minimum output entropy with a high precision. 



A. Iterative minimization of the output entropy 

We shall describe an algorithm that generates a se- 
quence of states \ipk), k = 0, 1, 2, . . ., such that 

• The output entropies Sk — S($(\ipk+i){'tpk+i\)) are 
non-increasing, i.e. Sk+i < <Sfc. 

• The algorithm converges to an extremal point of a 
functional S($(|V>)(t/>|)). 

The algorithm exploits the following simple observa- 
tion. Consider a functional 

fc(p,»j) = -Tr[*(p) log $(77)] 

that depends upon a pair of density operators p and 77. 
Let us show that the global minimum of h coincides with 
S m i n . Indeed, let us use an identity 

h( P , 77) =sm P ))+sm P )\ ISM), (42) 

where 5(u;||cr) = Tr[w(logw — logo - )] is the relative 
entropy. Since the relative entropy is non-negative, 
5(w||cr) > for any u>,a, and 5(w||ct) = iff u> = a, 
we conclude that 

S m in(®) — mm/i(p, 77). (43) 



The key observation is that finding the minimum of 
h(p,r]) with respect to each individual variable is an ex- 
actly solvable problem. Indeed, if one fixes 77 then 

min h(p, 77) = — maxTr [p <I>*(log <&(t7))] 
p p 

= -A mox ($*(log$(77)), (44) 

where $* is the linear map conjugated to $ (with respect 
to the trace inner product) and X max (X) is the maximum 
eigenvalue of an operator X. The optimal state p is ob- 
viously the highest eigenvector of $*(log<i>(?7)). On the 
other hand, if one fixes p and minimizes h(p, rj) over 77, 
the absolute minimum is achieved at p — rj, as immedi- 
ately follows from Eq. 1(4*21 . 

We shall attempt to find the minimum of /?,(p, 77) by 
earring out a sequence of alternating one- variable mini- 
mizations. This iterative minimization procedure gener- 
ates a sequence of pure states |t/j ), \ipi), ■ ■ ■ , such that (i) 
the initial state \ipo) is chosen randomly; (ii) \ipk+i) is the 
highest eigenvector of an operator <E>*(log $( (?/>/; |)) (if 
the highest eigenvalue is degenerate, \ipk+i) is chosen ran- 
domly from the corresponding eigenspace). Consider a 
sequence 

Sjt = S'($(|Vjt)(Vfc|)), fc = l,2,...,oo. 

Since each one-variable minimization does not increase 
the objective function, we have Sk+i < Sk, i-e., the se- 
quence {Sk} is monotone decreasing. 

Denote 5* = lim^oo Sk ■ Let be any limiting 
point of the sequence {IV^)} and p* = \ip*}{ip*\. Obvi- 
ously 

/i(p*,p*) < h(p, p*), and h(p*, p*) < h(p*, rj) 

for any states p and 77. Therefore a point p = 7/ = p* is an 
extremum of /i(p, 77). Thus p* is an extremal point of the 
output entropy 5($(p)) = h(p,p). Obviously, p* can not 
be a (local) maximum of h. However, it might be a local 
minimum or a saddle point. In principle, saddle points 
can be eliminated by computing a Hessian of h(p, p) at 
the point p*. 

It seems very unlikely that the iterations can be 
trapped in a very narrow and shallow local minimum, 
since one-variable minimization allows one to 'tunnel' 
through potential barriers. We can also expect that 
the method should work well for channels that are suffi- 
ciently noisy, since a landscape of the objective function 
S($(\il})(ip\)) becomes more flat as $ becomes less noisy. 
Symmetry of a channel also might be an important is- 
sue because one- variable iterative minimization typically 
works well for highly symmetric functions. 

Before applying the algorithm to fermionic product 
channels we have tested it on products of one-qubit uni- 
tal channels for which the exact value of S m i n is known, 
see (?). The algorithm always converges to the correct 
value of S m in for channels with sufficiently large amount 
of noise. For example, let O be a one-qubit Pauli depo- 
larizing channel, such that errors a x , a y , and a z occur 
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with a probability p each. Then the iterative algorithm 
correctly computes S min (0® 4 ) = 4H(1 — 2p) in a region 
p > 0.05. 

B. Numerical simulations 

In the rest of this section we discuss results of nu- 
merical simulations performed for several families of 
fermionic product channels on three and four qubits. 
For each particular channel the minimum output entropy 
has been evaluated using the iterative minimization al- 
gorithm with 64 iterations and 16 independent choices of 
the initial state. 

A plot on Figure 1 consists of two groups of points 
marked by '+' and 'x'. Each of these point corre- 
sponds to some particular 3-qubit channel from the fam- 
ily Eq. The points '+' correspond to a symmetric 
channel 

$+(c p ) = foc p , p = l,...,6, (45) 

where < b < 1 is a parameter specifying the channel. 
The points ' x ' correspond to a channel with a non-trivial 
distribution of the coefficients b p , namely 

$x(c p ) =bP/\ p =!,...,&. (46) 

The parameter b is plotted on the horizontal axis. On 
the vertical axis we plotted the smallest output entropy 
Smin found by the algorithm. The dotted lines represent 
the minimum output entropy that can be achieved on 
a Gaussian input state, i.e., the function S m in,g given 
by Eq. (|34|) . The inset plot shows the deviation S m in — 

Smin,g • 

It turns out that the algorithm always converges to 
the same value of S m in for the symmetric channel <&+ . A 
typical deviation S m in — S m in,g for this channel is about 
10 -9 which has the same order of magnitude as the nu- 
merical noise. 

As for the non-symmetric channel the algorithm 
always converges to the same value of S m i n unless b is 
close to 1. A typical deviation S m in — S m in,g is about 
10 -9 in the region < b < 0.7, so it is not shown on 
the plot. The deviation is always positive however. On 
the other hand, if b is close to 1, the iterations some- 
times converge to a local minimum with a value of S m i n 
exceeding S m in,g by several percents. 

Numerical simulations have been also carried out for 
analogous 4-qubit channels. A plot on Figure 2 shows 
the smallest output entropy found by the algorithm for 
a symmetric channel defined as 

$+(c p ) = 6c p , p=l,...,8, (47) 

and for a channel $ x with a non-trivial distribution of 
the coefficients b p , namely 

$x(c p ) =& p/4 , p=l,...,8. (48) 
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FIG. 1: Three-qubit channels. Vertical axis: the smallest 
output entropy Smin found by the iterative minimization al- 
gorithm for the channels and <& x - Horizontal axis: the 
parameter b specifying the channels, see Eqs. 145I46H . Dotted 
lines: the minimum output entropy achieved on a Gaussian 
input (Smin.g)- On the inset plot: deviation S m m ~ S m in, g 
for the channel $x • 
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FIG. 2: Four-qubit channels. Vertical axis: the smallest out- 
put entropy Smin found by the iterative minimization algo- 
rithm for the channels "I>+ and $ x - Horizontal axis: the pa- 
rameter b specifying the channels, see Eqs. I47I48B . Dotted 
lines: the minimum output entropy achieved on a Gaussian 
input (Smin, g )- On the inset plot: deviation Smin — S m in, g 
for the channel <I> X . 



Here < b < 1 is a parameter specifying the channel. 

The algorithm always converges to the same value of 
S m in for the symmetric channel <i> + . A typical deviation 
Smin — S m in,g for this channel is about 10 -9 . 

For the channel <I> X the algorithm always converges 
to the same value of S m in in a region b < 0.55. The 
deviation S m in — S m m,g drops down to 10 -9 in a region 
b < 0.6, so it is not shown on the plot. 

The numerical results seem to support the conjecture 
that the minimum output entropy is achieved on a Gaus- 
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sian input state, i.e., that S , m i„($) = S m j„ !9 ($) for all 
fermionic product channels. The fact that the algorithm 
always converges to the same value of S ra „ for the sym- 
metric channels $+ suggests that a landscape of the cor- 
responding objective function S(^+(\ip){ijj\)) is particu- 
larly simple. It might be an indication that S m in can be 
computed analytically for the symmetric channels. 
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